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Abstract 

Predicting liabitat suitability under climate change is vital to conserving biodiversity. However, current species distribution 
models rely on coarse scale climate data, whereas fine scale microclimate data may be necessary to assess habitat suitability 
and generate predictive models. Here, we evaluate disparities between temperature data at the coarse scale from weather 
stations versus fine-scale data measured in microhabitats required for a climate-sensitive mammal, the American pika 
{Ochotona princeps). We collected two years of temperature data in occupied talus habitats predicted to be suitable (high 
elevation) and unsuitable (low elevation) by the bioclimatic envelope approach. At low elevations, talus surface and 
interstitial microclimates drastically differed from ambient temperatures measured on-site and at a nearby weather station. 
Interstitial talus temperatures were frequently decoupled from high ambient temperatures, resulting in instantaneous 
disparities of over 30°C between these two measurements. Microhabitat temperatures were also highly heterogeneous, 
such that temperature measurements within the same patch of talus were not more correlated than measurements at 
distant patches. An experimental manipulation revealed that vegetation cover may cool the talus surface by up to 10°C 
during the summer, which may contribute to this spatial heterogeneity. Finally, low elevation microclimates were milder 
and less variable than typical alpine habitat, suggesting that, counter to species distribution model predictions, these 
seemingly unsuitable habitats may actually be better refugia for this species under climate change. These results highlight 
the importance of fine-scale microhabitat data in habitat assessments and underscore the notion that some critical refugia 
may be counterintuitive. 
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Introduction 

Antliropogenic climate cliange has already profoundly afTected 
range and community structure for many taxa [1-5]. Predicting 
future changes in species distributions and interactions is a 
challenge with great consequences for conserving biodiversity. 
However, realistic predictions and viable conservation plans 
depend on accurate forecasts of future habitat suitability. 

Species distribution models (SDMs) are the most common 
approach to predicting habitat suitability under climate change. In 
this approach, a species' bioclimatic envelope is statistically 
determined from its present distribution, and suitable habitats 
are then predicted where similar climatic conditions will occur in 
the future [6]. Correlative SDMs have been criticized for failing to 
incorporate important factors like dispersal, biotic interactions, 
adaptation and behavioral plasticity [7]. However, the predictive 
power of these models is improved when used in conjunction with 
other approaches such as phylogenetics [8], mechanistic heat 
transfer models [9] or eco-physiological parameters [10]. 

In spite of recent criticism, few SDMs account for fine scale 
microhabitat features, which can profoundly influence microcli- 
mates and therefore habitat suitability [6, 11, 12]. In fact, a recent 
meta-analysis demonstrated that average climate grid lengths in 



SDMs were 10,000-fold larger than the animals they study [13]. In 
addition, few studies exphcitly incorporate data from the periphery 
of a species' range, though these populations may be particularly 
informative of climatic tolerance or critical habitat features [14]. 

The concept of refugia from climate change has received a great 
deal of recent attention in the literature [15-18]. Refugia are 
"safe-haven" habitats where species can persist in times of 
environmental change and may offer hope for in situ persistence, 
particularly for species with poor dispersal capacity [15]. A key 
feature of refugia is that their microclimates must be relatively 
stable and must buflFer species against climate variability [17]. 
However, refugia can be diflTicult to identify, particularly when 
they are surrounded by areas of less suitable habitat [17]. Indeed, 
a key aim of conservation science is identifying and protecting 
habitat features that can serve as refugia through environmental 
changes. 

In this study, we examined how micro-refugia may allow a 
climate-sensitive species to persist in habitat that appears 
unsuitable by the bioclimatic envelope approach. American pikas 
(Ochotona princeps) are small mammalian herbivores that have 
become widely considered indicators of climate change [19]. We 
hypothesized that vegetation features may create favorable 
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microclimates at low elevations by insulating against summer heat 
stress, which is known to affect pikas. Uncovering microhabitat 
variation in seemingly unsuitable climates provides important 
information about the thermal tolerance, capacity for plasticity 
and potential vulnerability of a species to future climate change. 
This information will therefore be vital for refining SDMs and 
targeting critical habitat refugia for conservation action. 

Materials and Methods 



Study organism 

American pikas {Ochotona princeps; order Lagomorpha) are 
small mammalian herbivores that are typically distributed in high 
elevation mountains of western North America because they 
require short, cool summers and winters with extended snowpack 
[20]. Pikas are obligate talus specialists and are rarely found 
outside of rockslides and boulder fields that provide these suitable 
microclimates. In addition, pikas are extremely sensitive to high 
ambient temperatures [21,22] and possess a limited dispersal 
ability [23,24]. Unlike most alpine mammals, pikas do not 
hibernate during the winter but spend the short alpine summer 
amassing food caches called haypHes [25]. It has been hypothe- 
sized that warm summer temperatures may prevent pikas from 
constructing adequate haypiles, resulting in over-winter mortality 
[26]. This effect may be compounded by reduced snowpack, 
which exposes the animals to colder winter temperatures [27-29]. 

With climate change, parts of the pika's range are becoming 
unsuitable, resulting in upslope range retractions [27], population 
declines [28], and localized extinctions [27,30,31]. However, pikas 
persist in some regions with warm climates, including the 
Columbia River Gorge (CRG) in Oregon and Washington. Littie 
study has been devoted to CRG pikas, though they exist over a 
thousand meters below the species' previously recognized bioch- 
matic envelope [20,32]. Given this species' demonstrated sensitiv- 
ity to aspects of climate, marginal populations like the CRG have 
great potential to elucidate crucial features of micro-refugia and 
tractable conservation targets. 

Study sites 

We collected microhabitat data at four talus patches (ca. 15,000 
m~ each; hereafter, "sites") occupied by pikas in the CRG during 
June 2012-June 2014 (Table 1, Fig. SI). All sites were 32-35° in 
steepness, north to northwest facing and surrounded by a dense 
forest. One of the unique features of this region is a thick layer of 
moss that covers 25-80% of the surface of each rockslide. The 
moss also serves as a substrate for other plant cover, including 
graminoids, forbs and ferns. Sites ranged in elevation from 94 m 
to 437 m and varied in moss cover from over 65% ("high moss 
cover": sites 1 & 2) to less than 30% ("low moss cover": sites 3 & 4; 
Table 1, Fig. SI). Each of these sites is about a thousand meters 
lower in elevation than pikas are predicted to occur, based on the 
climate envelope approach [32]. 

We estimated potential solar exposure at each site using a 
previously-described insolation index [33]. The aspect and slope 
angle of each site were measured using a compass equipped with 
an inclinometer, and potential solar insolation was calculated as 
sine(slope) x cosine(aspect). This index ranges from —1 to 1, 
where values of 1 indicate steeper north-facing slopes with little 
solar exposure, and values of — 1 indicate steep south-facing slopes 
with high exposure [33]. 

We also sampled alpine microhabitats typical for pikas at two 
additional sites on the north face of Mt. Hood, approximately 
30 km from sites 1-4 (Table 1). Clast size was similar across CRG 
and Mt. Hood sites, with most rocks at each talus site having 
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Figure 1 . Pika microclimates compared to ambient temperatures. Temperature data were collected every two hours from June 201 2 to June 
2014 at four sites in the Columbia River Gorge. (A) Lines represent mean temperatures and shaded areas represent 95% confidence intervals. At high 
ambient temperatures, talus surface temperatures were lowest at sites with high moss cover. Forest surface temperatures were also cooler than 
ambient temperatures, but talus surface temperatures at sites with low moss cover closely tracked ambient temperatures. (B) Talus interstitial 
temperatures were functionally decoupled from variation in ambient temperature and remained a cool and constant 4-7 C at sites with high moss 
cover. (C) Time-trace of ambient, talus surface, forest surface, and talus interstitial temperatures measured at site 2 (high moss cover) during August 
2012. 

doi:1 0.1 371 /journal.pone.01 04648.g001 



dimensions in the range of 20-100 cm, which is preferred for pikas 
(Table 1). AH six sites were along publicly accessible trails on land 
owned by the United States Forest Service or Oregon Parks and 
Recreation Department. Appropriate permits for deploying 
dataloggers in Oregon State Parks were acquired from Oregon 
Parks and Recreation Department, permit no. 012-11. 

Macrocllmate temperature measurements 

We downloaded ambient temperature data of the type used for 
SDMs from the Western Regional Climate Center (available: 
http://www.raws.dri.edu/wraws/orF.html) at the Cascades Locks 
station (45°40'10"N, 12r52'54"W, elevation 128 m). This 
weather station is centrally located in the CRG, 10.6 km from 
our four low elevation field sites, on average. It is closest (3.4 km) 
to site 3. To further characterize ambient temperatures at CRG 
sites, HOBO Pendant temperature dataloggers (model UAOOl-08; 



Onset Computer, Bourne, MA) were suspended from a tree 
branch approximately 2 m above the ground and protected from 
solar radiation by a white plastic shield. All measurements labeled 
"ambient" reflect these 2 m shade air temperatures. 

Microclimate temperature measurements 

We measured temperatures in pika-relevant microhabitats with 
HOBO Pendant temperature dataloggers, which were housed in 
waterproof, plastic cases. At each site, we placed loggers at the 
talus surface in vegetated areas near the bottom of the talus slope 
where we had observed pikas foraging [34]. AH surface loggers 
were shielded from direct solar radiation by moss. To characterize 
sub-surface temperatures, we threaded a second, paired sensor 
into the talus interstices to reach a final depth (i.e., vertical distance 
directly below the surface logger) of 80-100 cm. In some cases, 
small sub-surface rocks were removed for sensor insertion and 



PLCS ONE I www.plosone.org 



3 



August 2014 | Volume 9 | Issue 8 | e104648 



Relevant Microclimates and Habitat Suitability 




1 10 100 1000 10000 

Distance (m) 

Figure 2. Spatial variation in pilca-relevant microclimates and 
ambient temperatures. Thermal semivariance is shown as a function 
of distance between dataloggers at the 4 CRG sites during June-August 
2013. Lines represent mean semivariance and shaded areas represent 
95% confidence intervals. Temperatures in pika-relevant microclimates 
are far more heterogeneous and less spatially correlated than ambient 
shade temperatures measured 2 m above the ground. 
doi:1 0.1 371 /journal.pone.01 04648.g002 

immediately replaced, but vegetation cover at the surface was not 
disturbed. Finally, at sites 2 and 3, we deployed a third logger in 
the surrounding forest to characterize a potential midday thermal 
refuge. For consistency, these forest surface loggers were also 
shielded from solar radiation by moss. 

Average summer temperature is known to impact pikas [27- 
29,33]. Furthermore, summer temperatures in the CRG are 
particularly unusual for this species [32]. To adequately capture 
average summer temperatures and to characterize spatial variation 
in summer temperatures, we deployed three additional surface- 
interstitial temperature logger pairs at each site during June- 
August 2013. To control for daily patterns of shade cover on the 
talus, one logger pair was placed near the top of the talus slope, 
and the other two logger pairs were placed approximately halfway 
up the talus on the east and west sides. Thus, at each site, we 
collected temperature measurements every two hours at four 
locations at the talus surface and four locations in the talus 
interstices during this time period. All datalogger locations were 
placed in areas of the talus where pikas were observed foraging, 
determined from behavioral observations [34]. 

Finally, to provide an initial characterization of humidity 
patterns between sites, we deployed HOBO relative humidity 
logger (model U12-011; one per site) at sites 2 and 3 during June 
2012-August 2013. Loggers were placed in plastic housing with 
two mesh walls that permitted sufficient air exchange for accurate 
humidity readings but prevented the logger from direct contact 
with substrates or precipitation. This housing was then threaded 
approximately 80 cm into the talus interstices. 

Spatial variation in temperature 

To characterize thermal variation at very small scales, we 
deployed additional surface, interstitial and ambient loggers at site 
2 and site 3 fromjuly 14— August 30, 2013 in close proximity (1 m) 
to the sensor network described above. We then characterized 



spatial variation in temperatures within and between sites with 
variograms, which are widely used in geostatistics to represent 
autocorrelation between measurements observed at different 
spatial locations. Variograms were recendy recommended as a 
spatially explicit tool for examining thermal variance in microcli- 
mates as a function of distance between measurements [13]. The 
semivariance (y) of temperature measurements at a given time 
point is half the average squared difference between logger-values 
(x), separated by a distance h, as given by: 

y(h)= — — (Xj — Xjf', where i, j are specific logger pairs 

separated by distance h, and N/^ is the number of logger pairs that 
are separated by this distance [35]. Lower semivariance therefore 
indicates higher autocorrelation between temperature measure- 
ments at a given separation distance. We computed semivariances 
with the variogram function in the R package gstat [36,37] for 
temperature measurements during June-August 2013. This 
function calculated semivariance at four average separation 
distances: 1 m (between closely placed loggers at sites 2 and 3), 
39 m (between loggers within a site), 3,622 m (between loggers at 
nearby sites) and 10,908 m (between loggers at distant sites). 

Vegetation cover and microclimate 

To provide an initial test of the hypothesis that vegetation 
insulates against extreme surface temperatures, we manipulated 
vegetation cover at an additional talus patch (ca. 1 km from site 1). 
A patch of moss approximately 1 m in diameter was experimen- 
tally removed and relocated nearby where moss did not naturally 
grow (Fig. S2). We then measured temperatures every two hours 
from May 28^une 13, 2012 at four locations: under an 
vmmanipulated patch of moss, under the transplanted moss, under 
a pile of bare control rocks that did not naturally have any moss, 
and under the pile of bare rocks where we removed moss (Fig. S2). 
The bare rocks in this experiment were approximately 15—20 cm 
in average dimension. All four locations were within 2 m of each 
other, and therefore should experience similar patterns of shade 
and sun exposure throughout the day. We compared daily average 
and maximum temperatures with two-way, repeated measures 
ANOVA, using treatment (i.e., manipulation vs. control) and 
cover (i.e., moss vs. rocks) as main effects and date as a repeated 
measure. 

Microclimate and elevation 

Temperature should decrease as elevation increases according 
to region-specific lapse rates, which may also vary by season [38]. 
In this region, previously reported lapse rates are — 2°C/m for 
minimum temperatures and — 7°C/m for maximum temperatures 
[39]. To investigate the impact of elevation on microclimate and 
calculate pika-relevant lapse rates, we placed paired surface- 
interstitial dataloggers at each high elevation site, which recorded 
temperatures at the same two-hour intervals as CRG loggers 
during June-August 2013. To test the effect of elevation on 
temperatures in pika-relevant microclimates, we computed 
Pearson correlation coefficients between elevation and average, 
maximum and minimum temperatures when all sites were free of 
snow. 

Data deposition 

The temperature and relative humidity measurements collected 
for this study are freely available in USPACE (http:/ /uspace.utah. 
edu), the University of Utah's institutional repository. The 
accession number for these data is ark:/87278/s6b583mp 
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Figure 3. The effect of vegetation cover on summer surface temperatures. (A) Average and maximum summer surface temperatures were 
significantly cooler under moss, compared to rocks. Both types of surface temperature were cooler than ambient temperatures. (B) Time-trace of 
temperatures during this experiment. 
doi:1 0.1 371 /journal.pone.01 04648.g003 



(available: http:/ / coiitent.lib.utah.edu/ cdm/ref/ collection/ uspace/ 
id/10610). 

Results 

Macroclimate temperature measurements 

We found no difference between 2 m shade ambient temper- 
atures collected at sites 1-3 and temperatures measured at the 
Locks weather station for daily averages {F}j28S — 0.654, p = NS), 
maxima (fjj2ss = 2.162, ^ = NS) or minima (f^j^ss = 0.654, 
p = 'NS) collected between October 2012 and August 2013. We 
therefore present data from our ambient loggers in future analyses 
because these data were collected at the same two-hour intervals as 
microclimate data. 

Microclimate temperature measurements 

We observed substantial divergence in microclimates relevant to 
pikas compared to 2 m shade ambient temperatures (Fig. 1). 
Although surface sensors were below moss cover at all sites, 



patterns in talus surface temperature were dependent upon overall 
moss cover at the site. At sites with high moss cover (>65% 
coverage), talus surface temperatures rarely exceeded 20°C, even 
at ambient temperatures over 35°C (Fig. lA, C). In contrast, at 
sites with low moss cover (£35% coverage), talus surface 
temperatures were only a few degrees cooler than ambient 
temperatures, on average (Fig. lA). Temperature at the forest 
floor was also 4— 9°C cooler at high ambient temperatures, but 
rarely dipped below 0°C on very cold days (Fig. lA). 

At all sites, interstitial temperatures were strongly divergent 
from ambient temperatures (Fig. IB). Even at ambient tempera- 
tures over 35°C, interstitial temperatures at sites with high moss 
cover remained a cool and constant 6-8°C, resulting in 
instantaneous disparities of up to 31.5°C between 2 m shade 
ambient temperatures (measured on-site) and microhabitats 
relevant to pikas (Fig. IC). Daily variation in interstitial temper- 
ature was < 1°C during the summer at sites with high moss cover, 
essentially decoupling interstitial temperatures from temperatures 
at the surface (Fig. IB, C). At sites with low moss cover, interstitial 



PLOS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | e104648 



Relevant Microclimates and Habitat Suitability 



o 

o 



30 



25 - 



-1— • 
CO 



20 - 



0 
Q. 

E 

.CD 



15 - 



10 - 



1., 



Talus Surface 
Talus Interstices 



1 I I I 

0 500 1000 1500 

Elevation (m) 

Figure 4. Average summer temperatures by elevation. Points 
represent mean summer temperatures at the talus surface (filled) and 
1 m into talus interstices (open). Bars range from mean daily maxima to 
mean daily minima in each of these habitats. To show bars that overlap, 
interstitial points are shifted right by 20 m at some sites. Dotted lines 
represent previously-calculated lapse rates for the Pacific Northwest, 
from [39]. 
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temperatures were more variable but remained 5-1 1°C cooler on 
average than the talus surface. 

Finally, relative humidity varied with ambient temperature (Fig. 
S3). At ambient temperatures below 10°C, both sites 2 and 3 had 
uniformly high relative humidity (>90%). At warmer ambient 
temperatures, relative humidity at both sites became more 
variable. However, site 2 (high moss cover) had consistendy 
higher relative humidity (75-95%) than site 3 (low moss cover; 55- 
75% humidity). 

Spatial variation in temperature 

Semivariance, or the average squared difference in temperature 
measurements separated by a certain distance, remained very low 
(<0.4°C^) among ambient temperatures, indicating high autocor- 
relation between these measurements at all scales (Fig. 2). In 
contrast, spatial variation in pika-relevant microclimates far 
exceeded spatial variation in ambient temperatures at all scales 
(Fig. 2). At a 1 m separation distance, semivariances for temper- 
atures recorded at the talus surface and in talus interstices were 
relatively low (<1°C^), indicating that temperature measurements 
separated by 1 m were highly correlated. However, at larger 
spatial scales, semivariance quickly increased to 12°C^ for 
temperatures recorded at the talus surface and 16°C^ for 
temperatures recorded in the talus interstices, indicating a very 
low degree of spatial autocorrelation. At separation distances 
larger than 1 m, microclimate temperatures recorded within the 
same site were not more similar to each other than temperatures 
recorded by sensors at the most distant sites. 

Vegetation cover and microclimate 

During the summer, average talus surface temperatures were 5- 
9°C cooler under both naturally occurring and transplanted moss, 
compared to bare rocks (Cover: F^/ = 239.2, j!)<0.0001; 



Fig. 3A). Surprisingly, there was also a significant effect of 
treatment, such that unmanipulated moss and rocks were also 
1-2°C warmer than manipulated moss and rocks (Treatment: 

F(i,i>) - 118.6, ^<0.001), but there was no interaction between 
cover and treatment. Maximum surface temperatures exhibited 
the same pattern (Cover: F(i }i) = 131.1, /><(). 0001; Treatment: 
■P'ri,/j} = 79.1, /)<0.001; Fig.'3A). Over the duration of tiiis 
experiment, both types of surface (moss and rocks) were cooler, 
on a\'(;rag(;, than ambient air temperature (Fig. 3A). Hcnvever, 
surface mi;asurements under rocks more closely trackc'd ambient 
temperatures; whereas surface measurements under moss were 
often up to 15°C cooler than ambient air temperature (Fig. 3B). 

IVlicroclimate and elevation 

Summer temperatures in the CRG were substantially cooler 
than predicted by region-specific lapse rates ([39]; Fig. 4). In fact, 
elevation was not predictive of talus interstitial temperatures or 
surface temperatures. All Pearson correlation coefficients between 
elevation and average, minimum or maximum temperatures were 
£0.6 and non-significant (p>0.2). The lowest elevation sites also 
had less than half the daily temperature range (i.e., daily 
maximum - daily minimum) of high elevation sites. 

Discussion 

We investigated the effects of vegetation and elevation on 
microclimates relevant to a sensitive species living in an atypical 
habitat. On the basis of climate envelope models alone, the 
Columbia River Gorge (CRG) appears to be unsuitable habitat for 
pikas because it does not have short, cool summers or long winters, 
which are commonly recognized as criteria for pika persistence 
[e.g., 20,32,40]. However, our results demonstrate that cool 
microclimates are available in the CRG and that these low 
elevation rockslides may actually be more suitable for pikas than 
the nearby alpine talus that is more typical for this species. 

Interstitial temperatures in the talus typically depend on heat- 
transfer processes at the surface (i.e., convection, conduction, solar 
and long- wave thermal radiation). However, during the summer in 
the CRG, daily fluctuations in this microhabitat are frequentiy less 
than 0.2°C and instantaneous disparities with ambient tempera- 
tures may be up to 31.5°C (Fig. IC), suggesting that, at some 
locations, the surface heat-transfer processes do not penetrate 
more than 1 m into the talus. As a result, interstitial microclimates 
are functionally decoupled from heat-transfer at the surface and 
are likely governed by other processes. CRG taluses have several 
features that are atypical for pika habitat, including a thick layer of 
vegetation covering the talus, high relative humidity, and reduced 
solar insolation (i.e., due to slope orientation and/or shade from 
the surrounding forest canopy). Each of these features may 
contribute to the significant disparities observed between the 
microclimate relevant to pikas and the macroclimate observed at 
weather stations. 

Vegetation features, and moss cover in particular, appear to 
buffer against climate variability, dHferentiating talus microcli- 
mates at some sites from the surrounding macroclimate (Fig. 1). 
Our high moss cover sites likely experience more shade (Table 1) 
and higher humidity (Fig. S3). These features may facilitate moss 
growth because mosses are predisposed to grow in shady, humid 
conditions [41]. However, our experimental results suggest that 
mosses may also actively buffer talus microclimates from ambient 
temperature extremes (Fig. 3). In other ecosystems and urban 
environments, vegetation can similarly lower surface temperatures 
of dark substrates by 5-20°C [42,43]. This cooling effect may be 
caused by changes in diffuse reflectivity (albedo), which can 



PLOS ONE I www.plosone.org 



6 



August 2014 I Volume 9 | Issue 8 | 6104648 



Relevant Microclimates and Habitat Suitability 



strongly afiect surface temperatures. Specifically, mosses may raise 
the albedo of the dark basalt that makes up the talus in this region, 
increasing the amount of solar radiation that is reflected. The 
albedo of basalt is 0.1, whereas albedo estimates of moss and 
lichen-dominated surfaces are two to four times higher [0.2-0.4; 
44]. For comparison, a change of albedo in 0.2 is capable of 
producing a 10°C change in midday surface temperatures [45]. In 
addition to increasing albedo, vegetation cover may also act as a 
layer of insulation [44], or it may have an evaporative cooling 
effect. 

Relative humidity and/ or precipitation may also play an 
important role in pika persistence and abundance in this habitat 
by maintaining high quality vegetation or water for evaporative 
cooling. During our study, the Locks weather station near our sites 
averaged ~150 mm of annual rainfall and 70-80% relative 
humidity. Similarly, relative humidity in the talus rarely dipped 
below 75% at our high moss cover sites (Fig. S3). Although 
differences in relative humidity may be responsible for the 
distribution of vegetation cover on talus in this region, the 
vegetation itself may also locally raise humidity through evapo- 
transpiration. Although we did not explicitly measure their effects 
on microclimates, sub-surface water or ice may also influence 
microclimate temperatures in this region. Water availability is a 
strong determinant of pika persistence and abundance in other 
parts of their range [30,46], and we regularly experienced cold air 
flowing out of talus cavities, which may indicate subsurface rock- 
ice or water features [47] . 

These low elevation sites also experience relatively litde solar 
radiation due to northerly slope orientations and daytime shade 
cover provided by the surrounding forest canopy. This reduced 
insolation likely has a large effect on both the magnitude and 
variability of temperatures in pika habitat. For example, although 
we did not directly measure solar radiation during this study, site 4 
likely received more sunlight because its surrounding forest canopy 
is shorter than other sites due to a wildfire in 1991 (Fig. SI). 
Indeed, although they are similar in elevation, slope angle, aspect 
and moss cover, summer surface and interstitial temperatures were 
approximately 5°C warmer and much more variable at site 4, 
compared to site 3 (Fig. 4). Similarly, a shorter forest canopy and 
higher insolation indices at our high elevation sites may have also 
contributed to greater temperature variability, compared to the 
CRG (Fig. 4). Unlike high elevations, CRG microclimates rarely 
experienced temperature extremes beyond the best estimates of 
thresholds for acute heat stress (25.5°C) [22] or cold stress (- 10°C) 
[27] in this species. 

Shade and vegetation cover are also likely responsible for the 
mild surface temperatures that we observed in the forest at sites 2 
and 3 (Fig. lA). Consistent with our findings, forest surface 
temperatures in other studies are typically 5-8°C lower than 
ambient temperatures measured above the ground [48,49]. Forest 
surface temperatures were typically a few degrees warmer than 
talus surface temperatures at sites with high moss cover (Fig. lA, 
C), but the forest may represent an important thermal refuge for 
pikas to remain active at midday at sites with low moss cover. 
Similarly, the forest may also serve as a winter refuge for CRG 
pikas, as forest temperatures also remained warmer than talus 
temperatures on the few ver^' cold days in 2012-2014 (<— 10°C; 
Fig. 1 A). Finally, thick understory vegetation in the forest may also 
make pikas less visible to predators or serve as an additional food 
resource [34]. 

Regardless of mechanism, the temperature disparities that we 
observed between weather stations and talus interstices in the 
CRG far exceed those in the microclimate literature. Previous 
studies examining ambient temperatures and microclimates at or 



below the soil surface report disparities up to 10°C between these 
measurements [12]. Similarly, burrow systems have a significant 
ability to buffer animals against temperature fluctuations, but 
instantaneous differences between ambient shade temperatures 
and internal burrow temperatures rarely exceed 10°C in the 
literature [50-55]. Rock-ice features in pika habitat in the Sierra 
Nevada range can also cool the average warm-season tempera- 
tures in talus interstices by 3.8°C, compared with the surface. In 
contrast, we observed diflerences of up to 11°C between average 
surface and interstitial microclimates in the warm season, and 
instantaneous disparities up to 31.5°C between talus interstitial 
temperatures and ambient temperature measured on site. For 
comparison, this disparity is roughly equivalent to 4500 m of 
elevation change (according to lapse rates for this region) and far 
exceeds the magnitude of observed and predicted directional 
climate change. 

CRG microhabitats were also highly heterogeneous. A lack of 
spatial autocorrelation between temperature measurements (even 
between loggers at the same site. Fig. 2) suggests that spatially 
uncorrelated noise exceeds deterministic trends in temperature at 
scales over ~10 m in these microhabitats. This result may be due 
to the patchy nature of moss and vegetation cover on the talus or 
to the movement of shade cover across the talus throughout the 
day. Interestingly, strong correlations were observed in the Great 
Basin between interstitial temperatures in montane pika habitat 
and ambient temperatures measured at valley weather stations 40- 
60 km away [29]. The temperature disparities and spatial 
heterogeneity that we report here highlight the uniqueness of 
habitats with high vegetation cover, and more generally, the 
refugial value of places where surface conditions can partially 
decouple ambient temperatures from sub-surface microclimates. 
Importanfly, high spatial heterogeneity in a habitat may also allow 
animals to behaviorally mitigate the impacts of directional climate 
change by shuttling between different areas within a habitat [12]. 

Taken t<)g("ther, these results suggest that tiu; low c-levation 
habitat in the CRG may actually be a better refuge from stresses 
caused by climate change than typical, high elevation habitat and 
that CRG pika populations have the potential to be relatively 
resilient to future environmental change. These results are a 
notable example of an unexpected or counter-intuitive thermal 
refuge, since the CRG was predicted to be unsuitable based on 
SDM results alone [e.g., 20,32,40]. Our results also suggest 
tractable conservation priorities for pikas and other thermally 
sensitive species in this region. Specifically, land managers could 
focus on protecting habitat features that contribute to the unique 
microclimates, including moss (e.g., from trampling by hikers or 
unsustainable harvest for the horticulture trade) and forest canopy 
cover (e.g., from logging or severe wildfires). 

An ecologically similar species that is likely to be direcdy 
affected by talus microclimates in this region is the bushy-tailed 
woodrat, Neoloma cinerea. Like pikas, these woodrats inhabit 
exclusively rocky habitats in montane forests and are sensitive to 
high ambient temperatures [56]. This species has experienced 
localized extinctions in response to past climate change [57] and 
appears to be suffering a range coUapse in response to 
contemporary climate change [1]. However, we have observed 
fresh woodrat sign at many talus patches in the CRG and 
speculate that cool microclimates and high vegetation availability 
also facilitate woodrat persistence at this unusually low elevation. 

Course-scale SDMs should not be discarded as a predictive tool, 
but climate data and their resolution must be selected with great 
care. We echo calls from Potter et al. [13] that the ideal spatial 
resolution for climate data in SDMs must be related to the body 
size of the species of interest. However, we also suggest that 
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heterogeneity of microclimate features and the spatial or temporal 
extent to which an organism samples these features should also 
inform selection of appropriate climate data in SDMs. For 
example, larger, less-mobile organisms are more likely to 
experience macroclimates observed at greater spatial grid lengths. 
Conversely, our results demonstrate that ambient temperature 
data collected at distant locations may have relatively Ktde 
relevance for small, mobile species that can rapidly shuttle 
between complex microhabitats, particularly in environments or 
seasons when surface features such as snow or vegetation decouple 
macro- and microclimates. Finally, the timing, frequency and 
amount of precipitation are rarely part of habitat suitability 
analyses, though they should be considered. 

Microhabitat data will be critical for developing and validating 
spatial statistics to understand the relationship between climate 
change and microclimate stability [13]. Relatively little is known 
about the degree to which microclimates wiU respond to 
macroclimatic changes or the time-scale on which these changes 
wiU occur. Collecting organism-relevant microhabitat data at finer 
spatial resolutions and longer time scales will be critical for 
bridging this gap and for d(;termining the resiliency of microcli- 
mate refugia to continued changes in broad-scale temperature and 
precipitation patterns. Such studies are becoming increasingly 
feasible, given the availability of relatively inexpensive and small 
temperature dataloggers with long battery lives. 

Our results highlight the need for fine-scale temperature 
measurements in microhabitats that are relevant to a focal species. 
When identifying refugia from climate change for conservation 
action, we can learn much from fine-scale observations of 
populations in seemingly marginal climates. Indeed, the conser- 
vation value of refugia lies in promoting species' ability to persist 
even under unfavorable climates [16]. 
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Supporting Information 

Figure SI Comparison of moss cover at four sites in the 
Columbia River Gorge. For reference, pikas are also shown at 
sites of high (E) and low moss cover (F). Photo credits: (A-E) J. 
Varner, (F)J. J. Horns. 

(TIFF) 

Figure S2 Moss transplant experiment. (A) In an area of 
the site where moss did not naturally grow, temperatures were 
measured under rocks and transplanted moss. (B) In an area of the 
site naturally covered with moss, temperatures were measured 
under rocks where moss was removed and a patch of 
unmanipulated moss. Arrows indicate datalogger locations. 
(TIFF) 

Figure S3 Relative humidity compared to ambient 
temperature. Relative humidity measurements were collected 
every 2 hours in the talus interstices at sites 2 and 3 during June 
2012-August 2013. Ambient temperatures (2 m height, shaded) 
were collected at site 3 at the same temperature inter\'als. Lines 
represent mean temperatures and shaded areas represent 95% 
confidence intervals. 
(TIFF) 
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